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Abstract 

We study the effect of lower dimensional geometry on the frequencies of the collective oscillations 
5h ' 

QJ of a Bose-Einstein condensate confined in a trap. To study the effect of two dimensional geometry 

: 

we consider a pancake-shaped condensate confined in a harmonic trap and employ various models 

,33 for the coupling constant depending on the thickness of the condensate relative to the the value 
P. 

'"j H of the scattering length. These models correspond to different scattering regimes ranging from 

'■ 

£h quasi-three dimensional to strictly two dimensional regimes. Using these models for the coupling 

o : 

, O p arameter and sum rule approach of the many-body response theory we derive analytical expres- 
^ H sions for the frequencies of the monopole and the quadrupole modes. We show that the frequencies 
of monopole mode of the collective oscillations are significantly altered by the reduced dimension- 

t> ; . 

<0 ality and also study the evolution of the frequencies as the system make transition from one regime 
G\ 

•O to another. 
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Recently, several theoretical and experimental studies devoted to the influence of dimen- 
sionality on the properties of a Bose-Einstein condensate in a harmonic trap have been 
reported in the literature QflBQBBBBE, 1Q, 11, 12, 13, 14, 15, lfi|. The reduction of 



dimensionality affects the physical properties of the condensates resulting in very different 
features from their three dimensional (3D) counterparts. In current experiments on Bose- 
Einstein condensates with alkali atoms confined in a magnetic trap the anisotropy parameter 
A (defined as the ratio between the frequencies of the trap in the z- and the transverse direc- 
tions) may be varied to achieve condensates with special properties which are characteristic 
of the low dimensionality. For example, by making the anisotropy parameter much larger 
than unity (A >> 1) a flatter and flatter (pancake-shaped) condensate can be produced. 
Such a pancake-shaped condensate is expected to exhibit special features of two dimensional 
(2D) condensate. It is well known that in the 2D case the scattering properties are very 
different as compared to the 3D case and this in turn lead to a significant modification of 
the boson-boson coupling constant. For example, the boson-boson coupling constant in 2D 
limit becomes density dependent even at a low density and zero temperature. In contrast to 
this for 3D case, to the lowest order in the density the interactions are described by constant 
coupling strength and deviations occur only when quantum depletion and finite temperature 
effects are taken into considerations jl7|. 

In a pancake shaped 3D condensate as the anisotropy parameter is increased the physical 
properties of the condensate first change due to the modified shape of the confinement 
12^ and then also due to the alteration of the scattering properties. With respect to the 
scattering properties of a flat pancake-shaped condensate three regimes can be identified 
which are characterized by different expressions for the coupling constant 10]. When the 
linear dimension of the condensate along z-axis given by a z = {K/muj z ) 1 / 2 (where uj z is the 
z-component of the angular frequency of the trapping potential and m is the mass of the 
trapped atoms) is still much larger than the 3D scattering length a (a z » a), the collisions 
still take place in three dimensions and this is referred to as quasi-3D (Q3D) regime. On 
further increasing the anisotropy parameter the condensate gets more tightly confined along 
the z-axis but the assumption that the scattering is unaffected begins to break down when 
a m a z and the condensate is said to be in quasi-2D (Q2D) regime. A fully 2D condensate 
is achieved when a z << a so that the collisions are restricted only in the transverse {x,y} 
plane. It is natural to expect that as the condensate evolves from a fully 3D to a fully 
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2D regime its static and dynamic properties undergo dramatic changes. In Ref. Q| the 
evolution of the density profile of a condensate has been investigated as the system crosses 
from the 3D to the 2D regime. They have found that the widths of the density distribution 
crucially depend on the collisional properties or the boson-boson interaction parameters in 
the different regimes. 

The main aim of this letter is to study the effect of lower dimensionality on the the 
frequencies of the collective oscillations of a condensate. In particular we calculate the 
collective oscillation frequencies of the quadrupole and monopole modes of a flat pancake- 
shaped condensate and study how these frequencies evolve as the system undergo transitions 



from the Q3D to the Q2D and to a strictly 2D regime. 
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or this purpose we make use of the 
By employing this method we 



sum rule approach of many-body response theory 
derive analytical expressions for the frequencies of the quadrupole and the monopole modes. 
In the following we first derive these expressions and then discuss the results. 

We consider a dilute condensate with N bosons confined in an anisotropic (pancake 
shaped) harmonic trap characterized by the frequencies u± and u z = Xu± with the 
anisotropy parameter A being much larger than unity. Within the density functional the- 
ory, the ground properties of a condensate can be completely described by the ground state 
condensate density p(r) in {x, y} plane. The ground state density of the condensate can be 
determined by minimizing the local density energy functional 

w_ 

2m 1 



E[p] = / d 2 r 



|V v / pW| 2 + ^(r)p(r) + e(p)p(r) 



(1) 



where v ex t(r) is the external harmonic potential in the transverse direction given by 



1 

Vext(r) = -mu 2 ± (x 2 + y 2 ) . (2) 

In the above equation (P) the first and the third terms represent the kinetic energy and 
the energy due to the interatomic interaction within local density approximation (LDA) 
respectively. Within this approximation the interaction energy per particle e(p) is given by 

e(p) = fp(r). (3) 

Where g is the coupling constant whose form depends on the collisional properties of the 
condensate. For example, g is independent of the density for the 3D case, on the other 
hand in the purely 2D regime it depends logarithimically on the density. In the following 



we briefly describe the models for the coupling constant g valid in the different collisional 
regimes. 

For the 3D system the coupling parameter g which is a constant completely determined 
by the s-wave scattering length a and it is given by 

9 = a. 4 

m 

When the linear dimension a z of the condensate along the z-direction is much larger the 
3D scattering length a (a z >> a), the collisions still take place in three dimensions. Under 
this condition the effective coupling constant (?Q3D which includes the effects of reduced 



dimensionality only is given by 



4. 
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g Q3D = 2^— (5) 
ma? 



On further increasing the anisotropy and and a z becoming comparable to a (a z « a), the 
collisions start getting affected by the tight confinement along the z-direction. Under such a 
condition the condensate is said to be in Q2D regime. The coupling constant in this regime 
is given by 

9Q2D l + 7 ^\ln{2{2^p{v)aa z )\ U 



Here p(r) is the ground state density of the condensate. The above expression was originally 
derived by Petrov et al. 0, 0] by studying the scattering properties of a bosonic system 
which is trapped harmonically trapped in the .direction and uniform in the {,,„} plane. 
A similar expression was later derived by Lee et al. |9| also by employing many-body T- 
matrix approach. It is important to note here that in the Q2D regime the coupling constant 
becomes dependent on the density in accordance with the behaviour of collisions in two 
dimensions. Finally, when a z becomes much smaller than a (a z << a), the collisions can 
be safely assumed to be taking place in two dimensions resulting in a 2D condensate. The 
coupling constant in 2D regime is given by 

g 2D = AnH2 1 ( 7 ) 

m \lnp(r)a 2 \ 

Notice that the expression for g2D also depends on the density p(r) but the information 
about the confinement direction is absent as it corresponds to the purely 2D case. The 
expression for g 2 n was derived in Ref. for a homogenous bose gas of hard disc. This 



4 



brm of g2D has been employed to study the properties cofined bosons in two dimensions 

aa 

and the rigorous justification for this use was provided by Lieb et al. [8j]. 
Having described the different models for the coupling constant, we now briefly discuss the 
method for obtaining the frequencies of the monopole and the quadrupole modes of collective 
oscillations. As mentioned earlier for this purpose we employ the sum-rule approach of the 
many-body response theory. The most important advantage of this method is that the 
calculation of frequencies requires the knowledge of the ground-state wave function (or the 
corresponding ground-state density) of many body system only. In accordance to the basic 
results of the sum-rule approach [3, I^J the upper bound of the lowest excitation energy is 
given by 

m ex = (8) 

V mi 

where 

m p 



^|(0|F|n)| 2 (^ n) P (9) 



is the p-th order moment of the excitation energy hu>o n associated with the excitation oper- 
ator F and Q ex is the frequency excitation. Here hjj n Q = E n — E is the excitation energy 
of eigenstate |n) of the Hamiltonian H of the system. The upper bound given by Eq. (JBJ) is 
close to the exact lowest excited state when this state is highly collective, that is, when the 
oscillator strength is almost exhausted by a single mode. This condition is satisfied by the 
trapped bosons in most of the cases. Moreover, Eq. (jHJ) can be used for computation of the 
excitation energies by exploiting the fact that the moments vn\ and m 3 can be expressed as 
expectation values of the commutators between F and H in the ground state |0) Q : 

mi = ~<0| [F\[H,F]} |0>, 

m 3 = i(0|[[Ft,F],[[F,[iJ,F]]]|0). (10) 



^ove equation we need to first 
lHj l the excitation opertaor F 



For the purpose of calculation of m\ and m 3 as given by the a . 
choose an appropriate excitation opertaor F. Following Ref. 
is written as 

F = x 2 + ay 2 (11) 

with a = 1 and a = — 1 for the monopole and the quadrupole modes respectively. By 
using the energy functional given by Eq. (P) along with the expression for e(p) we find after 
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some tedious although straightforward algebra following expressions for the frequencies of 
the quadrupole 

/ j.\l/2 

n q = V2[i + -) , (12) 



and the monopole 

^ = v^(l + ^ + ^) 1/2 , (13) 

modes. In the above equations T and U denote the kinetic and the harmonic confinement 
energies and they are given by 

U = ±yl J <i 2 r (x 2 + y 2 ) p(r) (14) 

On the other hand Yi nt arises from the interaction energy (third term) term of Eq. ([Tjl. 
As a result of this the values of Y int depend on the model of the coupling constant. The 
expressions for Y int (in the unit of Nhuj±) corresponding to the three regimes are given by 



yQ3D _ rpQ3D 

int int 

Y£ = (« + 2/r-2/ 2 2D ) (15) 



with k = a/ y/ft and the dimensionless interaction parameter a = a/a z . The general expres- 
sion for the interaction energy E\ nt {% = Q3D, Q2D, 2D) is given by 

EL = j rf 2 r| P 2 (r) (16) 

By using coupling constant given by Eqs. (JSJ)- ©in the above integral we obtain the inter- 
action energies in the corresponding regime. On the other hand, the expressions for other 
terms appearing in Eq. (JT5|l can be written as 



l + k\ln (2(2vr) 3 / 2 iVp(r)f) | 2 

jQ2D = ^ m f rfr 2 p 2 (r) 

2 "~ J 1 + k\ln (2(27r) 3 / 2 iVp(r)f ) | 3 

I™ = 2,N f dTV ^ 



\ln (Np(r)a 2 ) 



IT = 2nN I - ,t 2 ftll3 ( 17 ) 
' \ln(Np(r)a 2 ) | 3 



The density appearing in the above equation (Eq. ([17])) are normalized to unity and notice 
that the anisotropy parameter A is explicitly appearing in the integrals for the Q2D case. 
Before proceeding further we note that the expressions give n by Eqs. (|12j) and (|13|) for the 
Q3D case match with the corresponding results of Ref. 13[. 

It is evident from Eq. (fT2*|) that the frequency of the qudrupole mode is not explicitly 
dependent on the interaction energy and it is true in all the three regimes considered in 
this paper . We wish to note here that the Eq. (|T2*|) is identical to the expression for the 
requency of the quadrupole mode of a 3D condensate with constant coupling parameter 
In the absence of inter particle interactions we have T = U and this leads to harmonic 
oscillator result Q q = 2. On the other hand, when interaction energy is much larger than the 
kinetic energy so that the kinetic energy can be neglected ( Thomas-Fermi approximation) 
we get il q = y/2. In the general case one needs to know the value of kinetic energy of the 
ground state for accurate estimation of the frequency of quadrupole mode. 

Now we focus our attention on the frequency of the monopole mode of the collective 
oscillations. In contrast to the quadrupole case the frequencies of monopole mode are ex- 
plicitly dependent on the interaction energy. Therefore, the frequencies of the monopole 
mode will be affected by different models of the coupling constant. To illustrate the de- 
pendence of the interaction energy on the frequencies more clearly we substitute the virial 
relations associated with the three models 



T-U + Effi? = 



?Q3D 
J int 

T-U + E^ t D + M? 2D = 

T-U + E^ + I™ = (18) 

in the respective expressions in Eq. (|TB*|) to obtain 

t T Q2D T Q2D \ X l 2 

( T Q2D r Q2D\ 1 / 

^ = >{ 1 + -w--ir) (19) 

The first of the above equation shows that unlike the Q2D and the 2D cases the frequency 
of the monopole mode in the Q3D regime is independent of the coupling constant. The 



frequency of the monopole mode in the Q3D regime is given by the frequency of a 2D 
harmonic ocsillator. It is important to point out that in contrast to the Q3D case the 
frequency of the monopole mode of a 3D condensate is explicitly dependent on the interaction 
strength j2l|. 

Now we turn to the detail study of the frequencies of the monopole mode for the Q2D 
and the 2D models of the coupling constant. To this end we first need to evaluate the 
integrals given in the Eq. (JT7jl . For this purpose we employ the ground state densities 
p(r) of the condensates which are obtained within the Thomas-Fermi (TF) approximation. 
Furthermore, to obtain the densities within the TF approximation the spatial dependence 
of the coupling constants is also neglected by using results of the homogenous s yste m to 
relate the density to the chemical potential. It has been shown in the Refs. P.llo| that 
for large N the TF approximation and spatially independent form of the coupling constant 
yield sufficiently accurate results. 

We begin the discussion of the results with the values of relevant parameters from the 

n 

experiments of Gorlitz et al. with 23 Na atoms. These parameters are N = 10 5 , A = 
26.33 and a = 3.8 x 1CT 3 . We note here that this value of N is consistent with the TF 
approximation. These parameters indicate that the condensate produced in the experiment 
falls within the Q3D regime. The numbers obtained by us for the monopole frequencies of 
this system are VLq 2 d = 2.001 and VL 2 d = 2.084. These results clearly show that for the 
above system the Q2D result is very close to the corresponding Q3D number. However, the 
2D model overestimates the frequency of the monopole mode and this is anticipated as this 
model is not applicable to the condensate achieved in the above-mentioned experiment. 

Now to study the effect of different models of the coupling constant and their applicability 
we choose three different values of the parameter a: a = 3.8 x 10~ 3 , 0.33 and 2.68. These 
values are chosen such a way that the first, the second and the third numbers fall in the Q3D, 
the Q2D and the 2D regimes respectively. In addition to this we also choose A = 2 x 10 5 so 
that the condensate has negligible length in the z-direction and the motion along this axis 
is completely frozen. The results with these parameters are presented in Table I. Again we 
can see from Table I that for a = 3.8 x 10 -3 the numbers predicted by the Q3D and the 
Q2D models are very close and the corresponding result from the 2D model is quite higher 
than these two numbers. For a = 0.33 the numbers obtained by the both Q2D and 2D 
models are markedly different from the result of Q3D model. As has been discussed before 
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for a = 0.33 along with the large value of A the scattering properties start to get influenced 
by the confinement. In this situation it is expected that the Q2D and 2D models will give 
significantly different numbers in comparison to the corresponding Q3D result. In the light 
of our earlier discussion we expect that for a = 0.33 the coupling constant is better described 
by the Q2D model. On the other hand, for this value of a the frequency of the monopole 
mode obtained with the 2D model is higher than that of Q2D model. As the the interaction 
parameter a is further increased to a value a = 2.68 the scattering properties become truly 
two dimensional, consequently for this value of a the 2D model should be able to predict 
the frequency of the monopole mode of a two dimensional condensate. In contrast to the 
case of a = 0.33, for a = 2.68 the value of the frequency obtained by the 2D model is lower 
than the corresponding number from the Q2D model. 

Finally, for the sake of completeness we plot in Fig. 1 the frequencies of the monopole 
mode obtained with three different models as a function of the interaction parameter a. 
The curves are drawn with the anisotropy parameter A = 5 x 10 5 and the number of atoms 
N = 10 5 . It can be clearly seen from Fig. 1 that the frequencies of the monopole mode 
obtained with the three models of the coupling constant exhibit different trend with the 
increase in the interaction parameter a. For example, in contrast to the constant value of 
the frequency for the Q3D case the frequencies obtained by employing the Q2D and the 
2D models increase as a is increased. As mentioned before it is only for a « 1 the values 
of frequencies obtained by the Q3D and the Q2D models are very close and the 2D model 
gives quite different numbers. On the other hand, when a exceeds the value 10~ 2 , the effects 
of reduced dimensionality start affecting the scattering properties and both Q2D and 2D 
models give different results as compared to the Q3D numbers. Therefore, we conclude from 
our results that the collective frequencies of the monopole mode and their behaviour can 
be used to identify the dimensionality of the systems as the system evolves from Q3D to a 
strictly 2D regime. The results of theoretical calculations can also be used to test the validity 
of different models of the coupling constant by comparing them with the experimental values. 

In summary, we have calculated the frequencies of collective oscillations of the Bose- 
Einstein condensate confined in a flat pancake-shaped trap. The condensate is tightly 
trapped along z-axis such that the motion along this axis is frozen. For such a condensate 
depending on the value of the ratio between the size of the condensate along the z-direction 
and the s-wave scattering length, three different regimes can be identified identified. These 
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three different regimes are described by three different models for the boson-boson inter- 
actions. We have calculated the frequencies of the of collective oscillations corresponding 
to these three models. For this purpose we have used sum-rule approach of many body 
response theory along with the ground-state density obtained within the TF approximation. 
The main result of this paper is that the different models for the coupling constant are 
clearly manifested in the frequency of the monopole mode and they lead to distinct results 
in the region of their applicability. The effect of modification of the collision properties in 
two dimension start changing the monopole frequency when the value of a becomes more 
than 1CT 2 . On the other hand, quadrupole mode is not explicitly dependent on the coupling 
parameter and thus not affected by the boson-boson interactions. 
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Figure captions 



Fig. 1 Frequencies (in units of uj±) of the monopole mode of 5 x 10 5 23 Na atoms confined 
in a highly deformed trap with A = 2 x 10 5 as a function of the interaction parameter a. 
The solid line represents results for Q2D case while the corresponding 2D results are shown 
by the dashed line and the Q3D results are displayed by the horizontal line. 



12 



TABLE I: Frequencies of the monopole mode in the units of uj± for three different values of the 
dimensionless interaction parameter a calculated using Eq. H9t 'or N = 5 x 10 5 and A = 2 x 10 5 



a 




^Q2D 


^2D 


3.8x 10~ 3 


2.00 


2.001 


2.037 


0.33 


2.00 


2.049 


2.084 


2.68 


2.00 


2.292 


2.182 
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